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Abstract 
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Using the results obtained by the non commutative geometry techniques applied 
OO , to the Harper equation, we derive the areas distribution of random walks of length 

N on a two-dimensional square lattice for large N, taking into account finite size 
^ ■ contributions. 

Pacs Numbers : 05.45. +b, 72.15.Qm, 72.10.Bg 

O 

W> 

^jpjI Let us consider on a square lattice all closed paths of length N starting at the origin. 

For such a path T, let A(T) be its algebraic area. AsiVn oo, the average size of such a 
path increases as y/~N so that A(T) ~ N and the renormalized area will be a = A/N. We 
want to compute the probability distribution V(A, N) of the areas at large but finite N. 
In the limit N i— > oo, the distribution was computed first by Levy using brownian paths 
IjJ. We will give a method based upon the Harper model, allowing to compute the finite 
size corrections in a systematic way. The Harper model was designed in 1955 as the 
simplest non-trivial one describing the motion of an electron sitting on a two-dimensional 
square lattice and submitted to a uniform magnetic field. Let be the magnetic flux 
through the unit cell and let <fto = h/e be the flux quantum. We set 7 = 27r0/0 o - Then 
given m = (m 1 ,m 2 ) G Z 2 , we denote by W(m) the corresponding magnetic translations 
0. They satisfy the Weyl commutation rules 

W (m)W (m) = W(m + m!) exp(i—m A mf) , (1) 

where vnl A m = m' 1 m2 — m' 2 mi G Z. Note that 7 plays a role similar to the Planck 
constant in the canonical commutation relation. 
Harper's model is given by the following Hamiltonian : 

H=J2 W(a) , (2) 

|a|=l 
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where \a\ = \ai\ + \d2\ if a = (a±, a 2 ) £ Z 2 . In addition, one defines the trace per unit area 
as the unique linear map T on the algebra generated by the W(m)'s such that 

T(W(m)) = 5 mfi . (3) 

Then from (1), (2) and (3), we get : 



?(H N ) = E 



r : closed paths of length N 

where the sum is taken on the set of closed paths starting at the origin of length N. Note 
that N should be even to get a non zero sum. Let Qn be the number of such closed paths, 
we then get : 

^raax 

E Pn(A/N) exp(ixA/N) = fi^ 1 E exp(ixA(T)/N) = Q N l T{H N )\ 1=x/N . (4) 

From this relation we obtain : 

Q N = T(H N )\, =0 = J (2 cos k\ + 2 cos k 2 f 

aN+1 

= ^ (l + 0(l/iV 2 )) as N - 00 . (5) 

Moreover as N 1— > 00, for a given value of x, 7 = x/N tends to zero, so that we can use 
a semiclassical argument to compute T(H N ). It has been shown that the spectrum of H 
is made of Landau sublevels || : 

Etfr) = ± (4 - 7 (2£ + 1) + ^ [1 + {21 + l) 2 ] - 0( 7 3 )) £ = 0, 1, ... , 0(l/ 7 ) , (6) 

each with multiplicity per unit = 7/27T. So, 
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± e 



This gives : 



T(H m _ 4 " +1 */ 4 



27riV sinh(x/4) 
Using (4), we obtain the probability distribution : 

Vn{o) = J + 0{l/N) . (9) 
cosh (lira) 

We numerically computed Pat (a) from formula (8) and compared the result with exact 
numerical calculations (Figure 1). 
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FIG. 1. Scaling function for the probability of having a loop of area A for a random 
walk of N steps. Diamond, plus and triangle symbols correspond to the finite size data 
for N = 16, 18 and 20, respectively. Solid line corresponds to the universal function (9), 
as N i— > oo; whereas dashed lines include the 1/N correction term for N=20 and 40. 
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